1/5 


AD-A146  288 
UNCLASSIFIED 


DESIGNING  OPTIMAL  AOQL  SAMPLING  PLANS  A  COMPUTERIZED 
APPROACH(U)  FLORIDA  UNIV  GAINESVILLE  DEPT  OF  INDUSTRIAL 
AND  SYSTEMS  ENGINEERING  d  E  WALKER  ET  AL .  JAN  84 
RR-84-1  N00014-83-K-0093  F/G  12/1  NL 


III  10 

If  ia 

»-  m 

1 2.5 

III  ~~ 

|2-2 

Hill  1  1 

t  l> 

2.0 

llll  1  •  • 

= 

1.8 


MICROCOPY  RESOLUTION  TEST  CHART 


NATIONAL  K»*flAO  «'* 


AD- A 146  288 


!  l  l  ij . ;j»  -;  »  •  _  ^  1 

i<  I  i  d  1 

j  di-tsib'.iK.u  )  , : , .li i-.:?.  n  1 

• - - ' 

00 

CD 

o 

tH 

a 

© 

X 

DESIGNING  OPTIMAL  AOQL  SAMPLING  PLANS 
A  COMPUTERIZED  APPROACH 


Research  Report  Mo.  84-1 
by 

Jo  Ellen  Walker 
and 

Richard  S.  Leavenworth 


January  1984 


Department  of  Industrial  and  Systems  Engineering 
University  of  Florida 
Gainesville,  Florida  32611 


c 


,  ’V*. 

rM;v 


This  research  was  supported  by  the  U.S.  Department  of  the  Navy, 
Office  of  Naval  Research,  Under  Contract  N0014-83-K-0093. 

FINDINGS  OF  THIS  REPORT  ARE  NOT  TO  BE  CONSTRUED  AS  AN  OFFICIAL 
DEPARTMENT  OF  THE  NAVY  POSITION  UNLESS  SO  DESIGNATED  BY  OTHER 
AUTHORIZED  DOCUMENTS. 


ABSTRACT 


This  study  provides  an  algorithm  and  computer  program,  coded  in  Fortran 
IV,  to  develop  double  sampling  acceptance  plans  for  attributes  which  satisfy  a 
specified  Average  Outgoing  Quality  Limit  (AOQL)  and  minimize  the  Average 

■  :  ■  1  r  5  .  ’ , 

Fraction  -last-eeted -(-AFi)  at  a  specified  quality  level,  An  analysis  of  the 
steps  in  the  algorithm  is  provided  and  comparison  is  made  among  the  results  of 
program  runs  and  Dodge-Romig  AOQL  and  MIL-STD-1Q5D  plans. 


OBJECTIVE 


This  study  develops  a  computer  based  algorithm  to  derive  double  sampling 
plans  which  minimize  the  Average  Fraction  Inspected,  AFI(pl),  at  a  designated 
quality  level,  pi.  The  AFI  is  defined  as: 

AFI  =  1/N  ((L(p)*ns)  +  N*(l-L(p) ) :  Single  Sampling  (1) 

and, 

AFI  =  1/N  ( L(p)*nl)  +  Pa(n2)*n2  +  N*( l-L(p) ) ;  Double  Sampling  (2) 

L(p)  is  the  appropriate  formulation  for  the  probability  of  acceptance  as 
a  function  of  the  incoming  quality  level,  p.  This  study  employs  the 
binomial  distribution  which,  for  the  single  sampling  case,  yields: 

L(P)  -  l  (J)  Pd  (1-P)n*d  (3) 

d=Q  a 

where  n  is  the  sample  size  and  c  the  acceptance  number. 

For  Double  Sampling, 


1  n  di  n.-d. 
L(P)  =  l  ("1)P  1  (1-P)  1  1 

dx=0  al 


rr1 


c2"dl 

+  i  ,  i  i3;)i32)p 1  du-p) 

W1  d2=0  al  d2 


dl+d2, 


nl+n2"dl'd2 


(4) 


where  nl  is  the  first  sample  size  of  a  double  sampling  plan,  n2  is  the  second 
sample  size,  cl,  and  c2  are  the  acceptance  numbers  on  the  first  and  second 
samples,  respectively,  and  rl  is  the  rejection  number  of  the  first  sample. 

The  rejection  number  on  the  second  sample,  r2,  always  is  c2+l. 
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Pa(nl)  is  defined  as  the  probability  of  acceptance  on  nl  given  as 


1  n  di  n.-d. 

Pa(nl)  =  l  (.S1)  P  (!-P) 
d^O  1 


Pa(n2)  is  the  probability  of  acceptance  on  n2  given  as 


V  c2-dl 


Pa(n2)  =  l  L  Ajj 

dl=Cl+1  d2=0  1  2 


*  n  nn  di+d9  n  +n,-d.-d, 

i  Cotfjp1  2u-p)  1  2  1  2 


Feasible  sampling  plans  are  those  which  strictly  satisfy: 


L(pl)  >  1-a 


(5) 


(6) 


(7) 


where  pi  is  the  process  average  or  the  quality  level  considered  acceptable  for 
the  purpose  of  acceptance  sampling.  The  probability  of  acceptance  must  be 
least  1-a,  where  a  is  defined  as  the  Producer's  Risk,  (pi  and  a  are  specified 
design  parameter). 


2 


and  AOQL  <  AOQL*  (8) 

where  AOQL*:  specified  Average  Outgoing  Quality  Limit 

and  AOQL:  AOQ(p*)  =  p*( l-AFI(p*))  for  a  given  sampling  plan 

p*:  p  value  at  which  the  AOQ  function  reaches  its  maximum  value 


It  is  useful  at  this  point  to  differentiate  between  design  parameters  and 
acceptance  sampling  plan  parameters,  the  design  parameters  for  either  a 
single  sampling  or  double  sampling  plan  are  the  required  input  values,  pi 
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and  a,  the  so-called  Producer's  Risk  Point,  and  the  required  maximum  AOQL, 
AOQL*.  The  plan  parameters  are  the  outputs  of  the  programmed  algorithm,  ns 
and  c  for  a  single  sampling  plan  and  nl,  n2,  cl,  c2,  and  rl,  for  a  double 
sampling  plan. 
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PROBLEM  FORMULATION 

Initially,  a  single  sampling  plan  is  found  such  that  the  AFI(pl)  is 
minimized.  The  minimum  (ns,c)  combination  satisfying  AOQL  <  AOQL*  and 
satisfying  L(pl)  >  1-a  yields  the  optimal  single  sampling  plan  which  is 
designated  (n*,c*). 

The  bounds  on  the  parameters  of  the  double  sampling  plan  are  based  on 


these  results  from  deriving  the  single  sampling  plan. 

must  satisfy 

The  first  sample,  nl. 

c2  <  nl  <  n* 

(9) 

c2  >  c* 

(10) 

cl  <  C2-1 

(11) 

and  B  (nl,c2,pl)  >  1-a 

(12) 

where  cl  and  c2  are  the  acceptance  numbers  of  the  first  and  second  samples, 
respectively,  of  a  double  sampling  plan,  and  B(n,c,p)  is  the  cumulative 
probability  for  a  binomial  distribution  as  given  in  equation  (3).  For 
purposes  of  this  study,  rl  has  been  made  equal  to  r2,  that  is,  c2+l,  following 
the  Dodge-Romig  schemes.  The  second  sample,  n2,  must  satisfy: 

n2  >  n*-nl  (13) 

and 

BB(nl,n2,cl,c2,pl)  >  1-a  (14) 

where  BB(nl,n2,cl ,c2,p)  is  the  double  cumulative  probability  of  a  binomial 
distribution  as  given  in  equation  (4). 
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The  objective  is  to  minimize  AFI(pl)  subject  to: 

L(pl)  >  1-a 
AOQL  <  AOQL* 

For  each  feasible  cl,c2  combination,  an  optimal  nl,n2  combination  is 
such  that  the  AFI(pl)  is  minimized.  The  minimum  AF1  from  all  cl,c2 
combinations  yields  the  optimal  sampling  plan,  nl*,  n2*.  cl*.  c2*. 


found 
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ALGOROTHM 


The  algorithm  is  divided  into  two  parts,  one  for  single  sampling  and  one 
for  double  sampling.  The  bounds  on  some  parameters  of  the  double  sampling 
plan  are  based  on  the  results  of  deriving  a  corresponding  single  sampling 
plan,  initially  the  optimal  single  sampling  plan.  Thus,  to  simplify 
discussion  of  the  algorithm,  the  two  parts  are  treated  separately. 

Deriving  the  single  sampling  plan,  (ns,c). 

1.  Initialization:  ns=l,  c=0.  Input  a  (  >  0  ),  pi,  AOQL*  (  >  pi). 

2.  Compute  B(ns,c,pl). 

3.  Test:  Is  B(ns,c,pl)  >  1-a  ? 

(a)  If  so,  go  to  step  4. 

(b)  If  not  so,  increment  c,  set  ns=c+l,  and  go  1  <  ,)  2. 

(If  the  first  value  of  8(ns,c,pl)  is  not  >  1-a,  then  no  larger 
value  of  ns  can  satisfy  the  constraint.) 

4.  Compute  the  AOQL  for  ns,c. 

5.  Test:  Is  the  AOQL  <  AOQL*? 

(a)  If  so,  go  to  step  6. 

(b)  If  not  so,  increment  ns  and  go  to  step  4. 

(c  <  ns  <  N) 

6.  Compute  B(ns,c,pl). 

7.  Test:  Is  B(ns,c,pl)  >  1-a  ? 

(a)  If  so,  go  to  step  8. 

(b)  If  not  so,  increment  c,  set  ns=c+l,  and  go  to  step  4. 

8.  Compute  the  AFI(pl).  Program  outputs  ns  and  c  along  with  the  AF I  (pi)  and 
actual  AOQL  after  double  sampling  algorithm. 
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9.  Set  n*=ns,c*=c,  go  to  double  sampling  algorithm.  (Note  -  in  future  cycles 
step  8  is  not  employed.) 

Deriving  the  Double  Samplinn  Plan  (nl,n2,cl,c2) 

1.  Initialization:  Set  c2=c* ,  cl=0,  nl=c2+l,  and  n2=n*-nl;  go  to  step  4. 
(B(nl=c2+l,c2,pl)  >  1-a  is  always  satisfied.) 

2.  Compute  B(nl,c2,pl). 

3.  Test:  Is  B(nl,c2,pl)  >  1-a  ? 

(a)  If  so,  go  to  step  4. 

(b)  If  not  so,  increment  c2,  set  c=c2,  set  n=n*,  go  to  step  2  of 
single  sampling  algorithm.  (If  the  first  value  of  B(nl,c2,pl) 
is  not  >  1-a,  then  no  larger  value  of  nl  can  satisfy  the 
constraint. ) 

4.  Compute  B8(nl ,n2,cl ,c2,pl) 

5.  Test:  Is  BB(nl ,n2,cl,c2,pl)  >  1-a  ? 

(a)  If  so,  go  to  step  6 

(b)  If  not  so,  increment  cl  and  go  to  step  4.  (No  larger  value 
of  n2  will  satisfy  the  constraint.) 

6.  Compute  the  AOQL  for  the  plan  (nl ,n2,cl ,c2) 

7.  Test:  A0QL(nl,n2,cl,c2)  <  AOQL* 

(a)  If  so,  go  to  step  8 

(b)  If  not  so, increment  n2  and  go  to  step  6.  (AOQL  values  decrease 
with  increasing  values  of  n2.) 

8.  Compute  BB(nl,n2,cl,c2,pl) 

9  Test:  Is  BB(nl ,n2,cl ,c2,pl)  >  1-a? 

(a)  If  so,  go  to  step  10. 

(b)  If  not  so,  go  to  step  9a. 
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9a.  Test:  Does  nl=n*? 

(a)  If  so,  increment  cl,  set  nl=c2+l,  n2=n*-nl,  and  go  to  step  2. 

(b)  If  not  so,  increment  nl,  set  n2=n*-nl,  and  go  to  step  2. 

10.  Compute  AFI(nl,n2,cl,c2,pl). 

11.  Test:  is  AFI(current  plan)  <  AFI(previous  plan)? 

(a)  If  so,  store  AFI  (current  plan)  and  its  plan  parameters 
for  the  set  value  of  c2  and  go  to  step  12. 

(b)  If  not  so,  retain  AFI  (previous  plan)  and  its  plan  parameters 
for  the  set  value  of  c2  and  go  to  step  13. 

12.  Test:  Does  nl=n*? 

(a)  If  so,  increment  cl,  set  nl=c2+l,  n2=n*-nl,  and  go  to  step  2. 

(b)  If  not  so,  increment  nl,  set  n2=n*-nl,  and  go  to  step  2. 

13.  Test:  Does  cl=c2-l? 

(a)  If  so,  go  to  step  14. 

(b)  If  not  so,  increment  cl,  set  nl=c2+l,  n2=n*-nl,  and  go  to  step  2. 

14.  Test:  Is  AFI(nl ,n2,cl ,c2,pl)  <  AFI(nl,n2,cl,c2-l,pl)? 

(a)  If  so,  retain  AFI  (current  plan)  and  its  plan  parameters, 
increment  c2,  set  c=c2,  n=n*,  and  go  to  step  2  of  the  single 
sampling  algorithm. 

(b)  If  not  so,  previous  plan  AFI  and  its  parameters  (nl,n2,cl,c2)  is 
optimal . 
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ANALYSIS  OF  THE  ALGORITHM 

The  L(pl)  decreases  with  increasing  values  of  n.  This  result  reduces  the 
computational  time  greatly.  As  nl  increases,  B(nl  2, pi)  decreases  (Figure 
1).  Once  the  L(pl)  constraint  ’*  violated,  c2  is  incremented  rather  than 
computing  B(nl,c2,pl)  for  all  nl  <  n*  for  the  given  c2.  Also,  since  the 
double  cumulative  probabilities  decrease  with  increasing  values  of  n2,  the 
next  nl  value  is  generated  once  BB(nl,n2,cl,c2,pl)  <  1-a  instead  of  continuing 
the  computation  for  n2  <  N-nl  (Figure  2). 

Table  1  perhaps  best  simplifies  further  discussion.  The  tabular  results  are 
from  a  computer  run  wherein  the  specified  parameters  are  pl=0.03,  a=0.05, 

A0QL=0. 05,  and  N=500.  Each  cell  (block)  contains  the  (nl,n2)  pair  and 
corresponding  AFI(pl)  values  such  that  L(pl)  >  1-a  and  AOQL  <  AOQL*  for  each 
feasible  (cl ,c2)combination ,  i.e.  for  cl  <  c2.  Table  2  shows  the  optimal  plan 
within  each  cell  (cl,c2  combination),  the  column  minimum  (fixed  c2),  and  the 
optimal  plan  selected  from  the  column  minimums. 

By  studying  the  results  as  displayed  in  Table  1,  the  following  con¬ 
clusions  are  reached: 

1.  In  each  cell,  the  AFI(pl)  increases  with  increasing  n2  and  fixed  nl 
(Figure  3).  Therefore,  it  is  necessary  to  find  only  the  first  nl,n2 

combination  satisfying  the  constraints  rather  than  testing  all  feasible 
values  ofn2.  The  large  list  of  the  AFI  values  for  all  values  of  n2  and  a 
given  nl  is  not  included. 

2.  In  each  cell,  the  AFI  is  a  convex  function  of  each  minimum  nl,n2  combina¬ 
tion  (first  feasible  nl,n2  combination).  The  minimum  sampling  plan  for 
each  cell  is  one  associated  with  the  minimum  of  this  function  (Figure 
4.)  In  Table  1,  the  minimum  sampling  plan  of  the  first  cell  in  the  c2=2 
column  is  nl=ll,  n2=16. 
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BB(nl=2,n2,  cl=0,c2=2) 


TABLE  1 


Schematic  of  Program  Outputs,  Double  Sampling 


c2= 

2 

c2= 

3 

c2= 4 

c2=5 

i 

I 

n  i 

9 

10 

U 

12 

13 

14 

20 

n2 

23 

20 

16 

14 

13 

12 

16 

AFI(pl) 

.073063 

.072404 

.067187 

.067326 

.07038 

.073129 

.085011 

nl 

7 

8 

9 

10 

11 

12 

31 

n2 

54 

41 

35 

30 

26 

25 

6 

AFl(pl) 

.076090 

.062998 

.059519 

.056563 

.054611 

.057869 

.091863 

nl 

7 

8 

9 

10 
11 
12 

41 

n2 

67 

54 

47 

41 

37 

35 

6 

AFI(pl) 

.065837 
.057350 
.054957 
.05  2956 
.052887 
.054955 

.102421 

nl 

7 

8 

9 

10 

11 

12 

si 

n2 

80 

66 

58 

52 

47 

45 

6 

AFI(pl) 

.061384 

.055319 

.053936 

.053722 

.053912 

.0565S1 

.17795 
cel  1 

nl 

n2 

AF I  (pi) 

nl 

n2 

Afl(pl) 

nl 

n2 

AFl(pl) 

nl 

n2 

AFI(pl) 

18 

15 

.082664 

17 

38 

.071300 

16 

94 

.083959 

15 

110 

.077789 

19 

12 

.08409  7 

18 

30 

.068552 

17 

53 

.063678 

17 

67 

.059314 

20 

9 

.084268 

19 

24 

.066407 

18 

43 

.060595 

18 

56 

.057129 

21 

7 

.085656 

20 

21 

.067226 

19 

37 

.059908 

19 

49 

.056860 

22 

5 

.086196 

21 

17 

.066042 

20 

32 

.059564 

20 

48 

.057501 

23 

4 

.089133 

22 

15 

.067319 

21 

28 

.059740 

21 

39 

.057907 

26 

0 

.091797 

31 

6 

.087214 

41 

6 

.098  053 

47 

10 

.108607 

nl 

n2 

AFI(pl) 

nl 

n2 

AFl(pl) 

nl 

n2 

AFl(pl) 

26 

48 

.084916 

26 

70 

.080637 

26 

89 

.077190 

27 

23 

.080726 

27 

41 

.074589 

27 

57 

.07079  7 

28 

17 

.081398 

28 

32 

.073855 

28 

47 

.070319 

29 

13 

.082497 

29 

26 

.073923 

29 

40 

.070588 

30 

10 

.083788 

30 

22 

.074  750 

30 

34 

. 070989 

33 

4 

.088  216 

41 

6 

.095250 

51 

6 

.110683 

nl 

n2 

AFl(pl) 

nl 

n2 

AFl(pl) 

37 

27 

.089033 

37 

46 

.085970 

38 

19 

.090008 

38 

45 

.086241 

39 

14 

.091277 

39 

29 

.087349 

40 

11 

.093011 

40 

24 

.088489 

41 

8 

.094413 

41 

20 

.099759 

46 

1 

.103962 

51 

6 

.109071 

nl 

n2 

AFl(pl) 

47 

32 

.102570 

I 

48 

21 

.  103  790 

49 

15 

.105285 

I 

50 

12 

.107183 

L 

51 

9 

.108927 

13 


The  minimums  for  each  cell  are  given  in  Table  2. 

TABLE  2 


nl 

n2 

AFI(pl) 

nl 

n2 

AFI(pl) 

nl 

n2 

AFI(pl) 

nl 

n2 

AFI(pl) 

11 

16 

.067187 

11 

26 

.054611 

11 

37 

.052887 

10 

52 

.053722 

18 

15 

.082664 

21 

17 

.066042 

20 

32 

.059564 

19 

49 

.056860 

27 

23 

.080726 

28 

32 

.073855 

28 

47 

.070319 

37 

27 

.089033 

37 

46 

.085970 

47 

32 

.102570 

Column  Minimum 

nl=ll  cl=0 

n2=16  c2=2 

AFI=. 067187 


Column  Minimum 

nl=ll  cl=0 

n2=26  c2=3 

AFI=.054611 


Column  Minimum 

nl=ll  cl=0 

n2=37  c2=4 

AFI=. 052887 


Column  Minimum 

nl=10  cl=0 

n2=52  c2=5 

AFI=. 053722 


Minimum  double  sampling  plan:  cl=0,  c2=4,  nl=ll,  n2=37,  AFI(pl)=0. 052887 

A0QL=0. 49464 

Corresponding  single  sampling  plan:  c  =  2,  n  =  26,  AF I (p  1 ) =0. 091796,  A0QL=0. 049770 
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AFI(pl) 


3. 


Each  column  gives  cells  for  all  feasible  cl  values  for  a  given  c2.  The 
cell  API  minimums  are  also  convex  in  cl.  Therefore,  the  minimum  sampling 
plan  for  a  given  c2  is  the  one  corresponding  to  the  minimum  of  these  cell 
minimums  (Figure  5).  For  example,  the  minimum  AF1  of  column  2  is  0.054611 
and  the  corresponding  sampling  plan  is  cl=0,  c2=3,  nl=ll,  n2=26. 

(Note:  It  is  not  necessarily  the  case  that  the  minimum  always  will  occur 
at  the  lowest  feasible  value  of  cl.  For  example,  for  the  parameters 
a=. 023606,  pl=.05,  A0QL*=.05,  and  N=35Q,  the  minimum  sampling  plan  is 
cl= 1 ,  c2=8,  nl=21,  n2=62.  However,  the  minimum  feasible  value  of  cl  is 
0.) 

4.  The  minimums  of  the  columns  also  follow  this  convex  pattern.  The  column 
minimum  AFI(pl)  values  in  Table  1  corresponding  to  c2=2,3,4,5  are 
0.072404,  0.0544611,  0.052887,  0.053722,  respectively.  This  provides  the 
minimum  samling  plan  as  well  as  a  stopping  criterion  for  the  algorithm 
(Figure  6).  Thus  the  minimum  of  the  function,  0.052887,  corresponds  to 
the  minimum  sampling  plan,  cl=0,  c2=4.  nl=ll,  n2=37. 

Single  Sampling 

The  first  feasible  ns,c  combination  is  found  such  that  B(ns,c,pl)  >  1-a. 
Since  the  binomial  probabilites  decrease  as  ns  increases,  the  first  infeasible 
value  of  ns  indicates  that  all  probabilities  associated  with  larger  values  of 
ns  for  a  given  c  are  <  1-a.  Thus,  once  a  point  becomes  infeasible,  c  is 
incremented,  and  the  process  begins  again. 

For  the  first  feasible  ns  value,  the  A0QL  <  A0QL*  constraint  is 
checked.  This  is  accomplished  in  a  subroutine  which  utilizes  the  Golden 
Section  Method  [3]  to  find  the  value  of  p,p*,  which  maximizes  the  A0Q 
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function,  i.e.,  the  AOQL .  p*  may  fall  in  the  interval  (0,1).  The  calculation 
of  the  required  initial  left  bound  (PLFT)  and  the  initial  right  bound  ( PRT ) 
are  based  on  this  allowable  interval.  In  each  iteration  of  the  search,  the 
AOQ  values  associated  with  PLFT  and  PRT  are  computed.  If  AOQ(PLFT)  a  AOQ(PRT), 
PRT=PLFT.  If  not,  PLFT=PRT,  i.e.,  after  comparing  the  AOQ  values,  only  one  of 
the  p  values  is  new.  Thus  only  one  additional  computation  of  the  AOQ  is 
necessary.  When  the  interval  on  p  is  less  than  0.0001,  the  search  ends,  and 
for  the  necessary  accuracy,  it  is  assumed  that  p*=(PLFT+PRT)/2.  If  AOQL  > 
AOQL*,  the  AFI  is  computed  and  the  process  begins  again. 

If  AOQL  <  AOQL*,  a  bisection  method  is  used  to  find  the  first  ns  value 
that  satisfies  the  AOQ  constraint.  The  initial  bounds  for  the  bisection  are 
ns  and  N.  If  the  AOQL  value  corresponding  to  N  is  greater  than  AOQL*,  the 
bisection  is  omitted,  c  is  incremented,  and  the  process  begins  again.  This  is 
allowable  since  the  AOQL  values  decrease  with  increasing  values  of  ns  (Figure 
7).  Use  of  the  bisection  method  in  this  case  appears  to  be  very  efficient. 

It  requires  iog(2)N  iterations  rather  than  the  N-ns  possible  iterations  needed 
using  a  total  enumeration  method. 

If  the  associated  binomial  probabilities  are>  1-a,  the  AFI(pl)  is  compu¬ 
ted.  Otherwise,  c  is  incremented,  ns  is  set  to  the  new  c  value,  and  the 
process  begins  again. 

Double  Sampling 

The  double  sampling  procedure  is  very  similar.  The  size  of  the  first 
sample,  nl,  ranges  from  c2+l  to  n*.  If  B(nl,c2,pl)  <  1-a,  c2  is 
incremented.  Since  these  cumulative  binomial  probabilites  decrease  with 
increasing  nl,  further  enumeration  for  a  given  c2  is  unnecessary.  Whenever  c2 
is  incremented,  new  bounds  on  n*  must  be  computed  using  the  single  sampling 
algorithm.  Since  c2  corresponds  to  c*.  a  new  n*  satisfying  the  single 

L 
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sampling  criteria  must  be  found.  The  double  sampling  algorithm  continues  with 
these  new  values  of  c2=c*.  c2  <  nl  <  n*,  and  n2=n*-nl. 


If  B(nl,c2,pl)  >  1-a,  the  corresponding  feasible  n2  >  ns*-nl  is  found. 

If  BB(nl,n2,cl,c2,pl)  <  1-  a,  cl  is  incremented  since  the  double  cummulative 
probabilities  decrease  for  the  first  feasible  nl,n2  combinations  of  a  given  cl 
(Figure  8).  Otherwise,  a  bisection  technique  similar  to  that  used  for  the 
single  sampling  case  is  used  to  find  the  lowest  n2  value  satisfying 
AOQL  <  AOQL*.  If  the  cumulative  double  probability  is  greater  than  or  equal 
to  1-a,  the  AFI(pl)  is  computed.  If  not,  nl  is  incremented  and  the  process 
begins  again.  For  a  fixed  nl,  the  first  feasible  (nl,n2)  pair  is  the  minimum 
of  the  AFI  function. 

The  minimum  sampling  plan  is  one  such  that  the  AFI(pl)  is  minimized. 
First,  the  minimum  for  each  cl  value,  the  cell  minimum,  is  determined.  This 
consists  of  the  nl,n2  combination  which  yields  the  smallest  AFI.  As 
previously  stated,  the  function  is  convex,  however  there  is  some  variation 
since  n2  does  not  decrease  in  uniform  increments  as  nl  increases.  In  Table  1, 
when  cl«l,  c2=3,  the  AFI  values  decrease  to  0.066407,  increase  to  0.067226, 
decrease  to  0.066042,  and  finally  increase  to  0.067319  as  nl  was  increased. 

If  the  function  were  strictly  convex,  the  computation  for  each  cell  would  end 
upon  the  first  increase  in  the  value  of  AFI(pl).  To  account  for  this  discon¬ 
tinuity,  the  algorithm  continues  for  five  additional  increments  of  nl  to 
reasonably  insure  that  a  global  minimum  AFI(pl)  is  reached. 

Next,  the  minimum  for  all  feasible  cl  values,  the  column  minimum,  is 
found.  Due  to  the  convexity  of  the  AFI  values  corresponding  to  the  cl  values, 
a  column  minimum  is  found  after  the  first  cell  minimum  increases. 

The  last  step  is  to  find  the  minimum  AFI  of  all  columns.  The  algorithm 
terminates  after  the  minimum  AFI  corresponding  to  a  c2  value,  (column 


22 


minimum),  is  larger  than  the  preceeding  one.  Thus,  the  minimum  sampling  plan 
for  the  design  parameters  is  found. 


The  computer  code  was  written  in  Fortran  IV  and  run  on  the  POP  11/34 
digital  computer.  For  further  clarification  of  the  computer  code,  flowcharts 
of  the  main  program  and  the  subroutines  that  compute  p*  and  AOQL  are  included 
in  the  appendix  as  well  as  a  listing  of  the  code. 


COMPARATIVE  RESULTS 


The  computer  program  was  run  with  various  AOQL  values.  The  following 
table  characterizes  the  effectiveness  of  the  program.  Sampling  plans  derived 
from  a  study  that  minimizes  the  ASN(pl),  plans  derived  by  Dodge  and  Romig,  and 
MIL-STD-105D  plans  are  shown  for  comparison. 

The  values  of  nl  and  n2  given  in  the  Dodge  Romig  tables  [1]  are  based  on 
the  largest  lot  size  in  each  range.  The  values  of  cl  and  c2  correspond  to  the 
mean  lot  size  and  to  the  mean  value  of  the  process  average.  This  ensures  that 
the  AOQL  values  over  the  range  will  not  exceed  the  specified  value,  however 
this  only  gives  the  average  most  economical  plan  within  a  range  of  process 
averages  and  lot  sizes.  The  computed  sampling  plan  will  vary  as  the  lot  size  and 
design  process  average  change.  For  example,  for  A0QL=1%,  a=0. 03465, 
pl=0. 0065,  and  a  lot  size  of  301  (lower  bound  of  the  Dodge-Romig  range),  the 
computer  program  gives  a  sampling  plan  of  cl=0,  c2=2,  nl=40,  and  n2=84  with  an 
AFI(pl)=0.  2170.  At  N=400,  the  upper  bound  of  the  lot  size  range,  cl=0,  c2=2, 
nl=46,  n2=73,  an  AFI(pl)=0. 1863.  Dodge  and  Romig  give  the  optimal  sampling  plan 
as  cl=0,  c2=2,  nl=55,  n2=60. 
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TABLE  III 


Comparisons  of  Computed  Results  with  Other  Standard  Plans. 


PUm  ualnp  Plant  uilnf  Do*|»-Bo»i|  MIL-BTD-103D  N1L-8T1M030 

«lnt«u«  AFI  niniaut  ASN  plan*  Nornal-ll  Tl|M«nW“II 

erlltrton  criterion  (Hln  AFI) 


AOOL-.Ol 

Cl-0 

Cl-0 

Cl-0 

cl-0  | 

mam 

•*  -.0347 

c2-2 

e2-2 

c2-2 

c2-l 

8  -.0736 

nl-43 

nl-36 

nl-33 

rl-2 

pl -.0069 

n2-77 

n2-3B 

n2-60 

r2-2 

*■2—2 

p2-.  039 

AFI-.2O03 

N-30 1-400 

n  1-30 

nl-80 

N— 330 

AFI— • 2321 

n2-30 

n2-QO 

AFI— •  2376 

AFI— .  4421 

AOQC-.  02 

cl-0 

Cl-0 

cl-0 

cl-0 

cl-0 

■*  -.0301 

c2-3 

c2— 3 

c2— 3 

c2-3 

c2-i 

O  -.3027 

nl-26 

nl-33 

nl-33 

r  1-3 

rl-2 

pi-.  013 

n2-33 

n2— 32 

n2-33 

r2-4 

r2-2 

P2-. 067 

AFI— » 1463 

N-30 1-400 

n  1  -32 

nl-32 

N-390 

AFI—. 184322 

n2-32 

n2-32 

AFI—. 1387 

AFI—. 2883 

AOOL-.  023 

cl-t 

Cl-I 

Cl-1 

Cl-t 

cl-0 

**  -.0389 

c2-4 

c2— 6 

c2-6 

C2-4 

c2-3 

/  -.  0044 

nl-42 

n  1-30 

nl-30 

rl-4 

rl-3 

.1-.  023 

n2-76 

n2— BO 

n2-60 

r2— 3 

r2-4 

.2— .09 

AFI-.  143113 

N-40 1-300 

n  1  -32 

nl-32 

ft*  300 

AFI-. 1838 

n2-32 

n2-32 

AFI -.0724 

AFI—.  1630 

AOOL-.OS 

€1-0 

Cl-1 

Cl-t 

cl-2 

Cl-1 

u  -.0433 

c2-3 

c2— 4 

C2-4 

c2— 6 

c2— 4 

>  -.1062 

Nl-ll 

Nl-23 

Nl-16 

n-3 

rl-4 

pi—.  04 

n2-46 

«2-23 

n2-34 

r2-7 

112-3 

p2-.  163 

AFI-. 078748 

N-401-3OO 

nt-32 

nl-32 

N— 430 

AFJ-.  0700 

n2— 32 

n2— 32 

AFI—.  0704 

AFI—. 1803 

AOOC-.03 

c  1—1 

el-l 

cl-0 

cl-3 

cl-3 

<  -.0442 

c2— 11 

c2— 4 

c2-4 

c2— 12 

c2— 1 1 

8  -.0630 

nl-22 

nl-26 

1*1-17 

rI-7 

rl-7 

pi—. 04 

na-iu 

n2— 24 

n3-33 

r2-13 

r2-l2 

p2— .  IB 

AFX -.03324 

N- 100 1—2000 

nl-00 

nl-BO 

N— 1300 

AFI— . 0631 

n2— BO 

n2-0O 

AFI—.  0673 

AFI—. 0982 

AOOL-.03 

cl-1 

c  1-0 

c  1— 2 

Cl-1  1 

•*  -.0403 

c2-4 

c2-4 

c2— 6 

c2— 4 

U  -. 0987 

nl-23 

nl-16 

r  1-3 

rl-4 

PI -.04 

n2— 44 

n2-24 

n2-33 

r2-7 

r2— 3 

p2— .  17 

AFI—. 072434 

N-30 1-4 00 

nl-32 

nl-32 

N— 330 

AFI—. 123807 

n2— 32 

n2— 32 

AFI—. 113002 

AFI—. 2043 
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The  computer  results  for  a 

2%  AOQL  with  a=0. 0478,  pl=0.02  and  N=1500 

follow: 

cl 

c2 

nl 

n2 

AFI(pl) 

0 

7 

28 

128 

.095858 

1 

7 

54 

144 

.093983 

2 

7 

81 

116 

.101264 

0 

8 

28 

193 

.096234 

1 

8 

55 

166 

.092043 

2 

8 

81 

142 

.098126 

0 

9 

28 

217 

.097668 

1 

9 

55 

190 

.091320 

2 

9 

82 

163 

.096105 

0 

10 

28 

242 

.101074 

1 

10 

55 

214 

. 092049 

2 

10 

82 

187 

.095498 

Dodge  and  Horn ig 

give  the  optimal 

plan  of  cl«l. 

c2=8,  nl=80,  n2=160  with 

the  AFI(pl)  being  0.1444.  The  minimum  sampling  plan 

found  above  (cl*=l,  c2*= 

nl*=55,  n2*=190)  has 

a  lower  AFl(pl)  value  than  that  corresponding  to  cl=l. 

c2=8.  The  difference  is  small; 

however  this  may  account  for  some  of  the 

discrepancies  in  the 

results. 

The  sampling  plans  derived 

from 

large  lot  sizes 

do  not  correspond  well 

with  the  Oodge-Romig 

plans.  For 

the 

parameters  A0QL=0.05,  a=0.0442,  pl=0.04 

and  N=1500, 

AFI  based 

Oodge-Romig 

MIL-STD-1050 

MIL-STD-1050 

algorithm 

AOQL  plans 

Normal -II 

Tightened-II 

pi  ans 

cl=l 

cl=0 

cl=5 

c  1=  3 

C2=ll 

c2=4 

c2=12 

c  2- 1 1 

nl=22 

n  1= 1 7 

rl=9 

Rl=7 

n2= 1 11 

n2=33 

r2=13 

R2=12 

nl=80 

nl=80 

b2=80 

n2=80 
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The  c*  value  of  the  single  sampling  plan  based  on  the  AFI  criterion  is  5. 

Since  c2  is  initially  set  to  c*,  the  acceptance  numbers  found  by  the  AFI 
algorithm  and  those  of  the  Dodge-Romig  plans  would  never  be  the  same.  It 
appears  that  for  large  values  of  N,  the  minimum  AFI  plans  resulting  from  the 
described  program  more  closely  resemble  those  of  the  MIL-STD-105D  plans. 

The  results  using  the  minimum  AF I { p 1 )  criterion  are  very  similar  to  those 
using  a  minimum  ASN  criterion.  The  AFI  based  plans  and  the  Dodge-Romig  plans 
deal  with  rectifying  inspection  whereas  the  minimum  ASN  based  plans  do  not. 

In  Table  III,  note  that  the  total  number  sampled,  nl+n2,  is  very  close  in 
the  AFI  algorithm,  ASN  algorithm,  and  Dodge-Romig  plans.  The  individual 
values  of  nl  and  n2  vary  to  a  larger  degree  than  the  total.  The  nl  value  from 
the  ASN  algorithm  plan  is  greater  than  that  of  the  Dodge-Romig  plans  and  n2  is 
less  than  the  n2  value  of  the  Dodge-Romig  plans.  Also,  nl(AFI)  is  less  than 
nl(Dodge-Romig)  and  n2(AFI)  is  greater  than  n2( Dodge-Romig) . 

The  MIL-STD-105D  tables  specify  a  fixed  sample  for  a  given  lot  size.  In 
most  cases,  the  total  number  sampled  does  not  correspond  well  to  the  other 
plans.  MIL-STD  Level  II  plans  for  Normal  and  Tightened  inspection  are 
included  in  the  table.  It  is  not  clear  which  specifications  most  closely 
related  to  the  other  sampling  plans.  In  some  instances  the  Normal  inspection 
plans  resemble  the  other  plans,  and  in  other  cases  the  Tightened  inspection 
plans  are  more  similar. 

The  computer  based  technique  is  advantangeous  in  that  plans  corresponding 
to  specific  AOQL  values,  lot  sizes,  and  process  averages  can  be  found  directly 
rather  than  relying  on  ranges  of  lot  size  and  process  averages  which  can  yield 
only  approximate  optimal  plans. 
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APPFNHIX 

SINGLE  SAMPLING  FLOW  CHART 


Initialization 
c=0 
ns=  1 


Compute  B(ns«c#pl)  using 
Subroutine  Probsl 


•^^1  s  B(nSiC.pl)£t -o£  ? 
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SUUROUTINE  SEARCH 


Initialization 
A  1=0 
A2«l 

T=.  00000001 
R  =  .  5*  ( SORTS— 1 ) 
H=A2-Ai 

PLFT=A1+(R*R)*H 

PRT=A1h(R*H) 


Compute  AOQ (PLFT) 
and  AQQ(PRT) 


Is  AOQ<PLFT> 
<  AOQ (PRT)? 


A1=PLFT 

H=A2-PLFT 

PLFT=PRT 

AOQ <  PLFT ) =AOQ  (  PRT ) 
PRT=A1+R*H 


!s  PRT— PI.  FT  <  T2. 


A2=PRT 
H=PRT— A1 
P RT=PLFT 

AOQ ( PRT ) =AOQ <  PLFT ) 
PLFT=A1+R*R*H 


:  s  PLFT— PLFT  <  TT. 


1 

P  r.  t«ir=(  PR  T+PLF  T)/P 
Compute  AOQ (Pc tar) 

-  •  -i 

Pstar=(PRT-t-PLFT>/2 
Compute  AOQ(Pstar) 

RETURN 

Compute 

AOQ (PRT) 

RETURN 

Compute 

AOQ (PLFT) 

QAFI . FOR 


0001 
0002 
0003 
0004 
0000 
0006 
0007 
0003 
0009 
OO 1  0 
001  1 
0012 
0013 
0014 
00 1  5 
0016 
00 1  7 
0013 
001  7 
0020 
0021 
0022 
0023 
0024 
0020 
0026 
0027 
0023 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0033 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
004  7 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 


C 

C*  tt*-**-***  ««  it*  »****!•  1- ***•«■*«•*  -tt**  **«*««»**«»### 

C  QUALITY  CONI  Li  ,  PROGRAM  TO  DERIVE  DOUBLE  SAMPLING 

C  PLANS  TO  MINIMIZE  AVERAGE  FRACTION  INSPECTED. 

C  : 

C  PROGRAMMED  H Y  >1.1  ELLEN  WALKER 

C  DEPARTMENT  OF  INDUSTRIAL  AND  SYSTEMS  ENGINEERING 

C  UNI  VERS  J  I  Y  01  FLORIDA 

C  GA1NESVI!  M  .  El.  OR  I  DA  32611 

C #  tt  turn-##  H-  S  »-«  -#■*  f 

r 

SUltRQUT  1  I'll  GAM 

INTEGER  CJ , C  3  <  C, CSTAR, C2M1 

INTEGER  K1.RJM1 

DOUBLE  PRECISION  SUMLOG 

BYTE  OUT  LI  I  CS) 

COMMON/BI  K4 /ALPHA. BETA 
COMMON/BLK6/C1, C2 
COMMON/Bt  K7 /SUMLOG  C  4000 ) 

COMMON/Bl K8/N 
C0MM0N/I3L  K12/0UTFIL 

C •»  «■•»■**•«■* ***»»**■&**  #*#■**  #■»■**•&■»•#■&*■«■*•«■•#■■***■«•■«•■«■*■*#■«■*■«■•«■■«■*•*  *■)(•■#■ 

C  INPUT  PARAMETERS 

C*  ** ******•«•*•«■•»  **•*  » 

KINDEX- 1 
WRITE  (5,  10) 

10  FORMAT! 1 OX,  'ENTER  VALUE  OF  ALPHA') 

READ  (5.  MAI  PHA 
WRITE  (5,  15) 

15  FORMAT! 1 OX,  'ENTER  VALUE  OF  PO  '  ) 

READ (5,  *)P0 
WRITE! 5, 20) 

20  FORMAT! 10X,  'ENTER  AOQL  VALUE') 

READ!  5,  *  )  AOCil 
WRITE! 5, 25) 

25  FORMAT! 1  OX,  'ENTER  LOT  SIZE') 

READ! 5, *)NNN 
C 
C 
C 

WRITE! 1, 23) 

28  FORMAT (' 1  ',/// 10X,  'DEPT.  OF  ISF  ' 

%/>  10X,  'UNIVERSITY  OF  FLORIDA  '/ 

*/SX,  5!'*')/  'DOUBLE  SAMPLING  SYSTEM  TO  MINIMIZE  AFI  ' 

*5! '*' ), 2X, /) 

WRITE (5, 30) ALPHA 
WRITE! 1, 30) ALPHA 

30  FORMAT! 10X,  'ALPHA-  ' ,  2X,  F8. 6) 

WRITE! 5, 35) PO 
WRITE ( 1, 35 ) PO 

35  FORMAT ( 1  OX ,  ' PO-  ' ,  2X «  F8. 6 ) 

WRITE!  5,  40  )  AOQL 
WRITE!  1,  40)A0QI 

40  FORMAT! I OX, ' AOQL - ' , 2X , F8. 6 ) 

WRITE (5,  4  5 ) NNN 
WR I TE ( 1 , 4  5 ) NNN 

45  FORMAT!  lOX,  'N  =  ',2X,I6> 
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OAF  I 


0058 
0059 
0060 
0061 
0061? 
0063 
0064 
0065 
0066 
006/ 
0060 
0069 
0070 
0071 
0072 
0073 
0074 
0073 
0076 
0077 
0070 
0079 
0000 
0001 
0032 
0003 
0084 
0085 
0086 
0087 
0088 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
0096 
0097 
0098 
0099 
OlOO 
0101 
0102 
0103 
0104 
0105 
0106 
0107 
0108 
0109 
0110 
0111 
01  12 
0113 
0114 


C*****************  ft************************************ 

C  COMPUTE  i>  •  MGLE  SAMPLING  PLAN 

C**************  if*************************************** 

c 

C  INITIALIZATION 

C 

C  #«■  ***•»■**»*■•&  *■  ■**■**■  ■#■•*■■#■*■**•*■*■*  ■*  **#  *  ******  •*••*#•  «•****  * 

NS-  1 

C^=0 

N=0 

AFI=1.  DO 
C 

WRITE! 5, 30) 

WR I TE  < 1 , 30 ) 

50  FORMAT! ///10X,  '***  SINGLE  SAMPLING  PLAN  *«•*■') 

C  *******  **■**■#■*■**** -a**  ■fr********************-**********  *** 

C  FIND  NS, C  COMBO  THAT  SATISFIES  L!PO>  G. T.  1-ALPHA 

C  ****************■*■******•#•**********■*■*•*■»■*•*■■*•**********■*'*■* 

55  CALL  PR0BS1 INS, PO, C, BXLEC) 

IF ! BXLEC.  LT.  !  1.  DO-ALPHA)  )C=C  H 
IF (BXLEC. LT. <1. DO-ALPHA ) ) NS-C+1 
C  ************************  ***  ***  ***  ****  ***  *****  ** ****** * 

C  SEARCH  TO  FIND  MIN  NS  VALUE  SUCH  THAT  ADQL  L.  T.  AOQL* 

C  ********  ***#**#**#»*##**»*****#MH«***»*#»#*******»*** 

60  CALL  SEARCH ( NNN, C, NS,  SADQ) 

IF ! SAOQ. LE. AOQL) GOTO  75 
NSTEMP=NNN 

CALL  SEARCH ( NNN, C, NSTEMP, SAOQ) 

IF ISAOQ. GT. AOQL ) C=C+1 
C 
C 

IF (SAOQ. GT. AOQL )NS=C+1 
C 
C 

BL=NS 
BH— NNN 

65  NSTEMP=I I  DINT ! ( BL+EH ) /2. DO ) 

CALL  SEARCH ! NNN, C, NSTEMP, SAOQ) 

IF ! SAOQ. LE. AOQL ) BH=NSTEMP 
C 

IF ! SAOQ. LE. AOQL ) AOQ=SAOQ 
C 

IF (SAOQ.  GT.  AOQL ) BL=NSTEMP 
IF< (BH-BL). EQ. 1. DO) GOTO  70 
GOTO  65 
70  NS=NSTEMP 

IF ( SADQ. GT. AOQL ) NS=BH 

C  ********** ******************************************* * 

C  CHECK  THAT  NS, C  COMBO  STILL  SATISFIES  L!P0)  CONSTRAINT 

C  **************  ***********  ***************************** 

75  CALL  PR0BS1 (NS, PO, C, BXLEC ) 

IFtBXLEC. GE. ( 1. DO-ALPHA) )GOTO  80 
C=C+1 
C 

NS=C+1 

C 

GOTO  60 
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0115  C  *****************  **  *******  *****************  ****  ******* 

0116  C  COMPUTE  API 

0117  C****** ******  a-**#**** ******  * ****************** *****  **** 

0110  80  ATIPO=NS*BXI  EC  t-NNN*!  i  DO-BXLEC) 

0119  AF  1  PQ=A'f  I PO/NNN 

0120  C 

0121  NS1ARNS 

0122  C. 

0123  WR  I TE  !  5<  05  )  NS 

0124  WR I TE ( 1 . 85 ) NO 

0125  85  FORMAT  (  /  /  1  OX .  'NS:-  '  ■  13) 

0126  WRITE <  5»  90 )C 

0127  WR  J  TE ( 1 j  90 )C 

0128  90  FORMAT  (  1  OX ,  'C-',I2> 

0129  WR ITE ! 5,  95 ) AF IPO 

0130  WR ITE! 1 (  95  > AF I PO 

0131  95  FORMAT ( J  OX ,  ' AF I ( PO > =  ' .  FO. 6 ) 

0132  C 

0133  100  CS  f  AR-C 

0134  C ************************************* ***************** 

0135  C  DOUBLE  SAMPLING 

0136  C  ****************************************************** 

0137  C 

0138  WRITE ( 5»  105) 

0139  WRITE! 1.  105) 

0140  105  FORMAT! ///10X,  ' **  DOUBLF  SAMPLING  PLANS  ** ' ) 

0141  PLAN= 1 . DO 

0142  C2“CST  AR 

0143  1 10  R 1=C2+1 

0144  C2M1=C2~1 

0145  ClPO=Cl+l 

0146  R 1M1— R 1—1 

0147  DDATI=NNN 

0148  TTCMIN=1. 

0149  C 

0150  C 1=0 

0151  JJ=0 

0152  C 

0153  WRITE ! 5.  115) 

0154  WRITE! 1.  115) 

0155  115  FORMAT! ///10X,  'C1',6X,  'C2',7X,  'Nl'.BX.  'N259X,  '  AF  I  ' ,  /  /  ) 

0156  C****************************************************** 

0157  C  CALCULATE  FIRST  SAMPLE  NUMBER 

01 58  C 

0159  C  FROM  RESULTS  OF  PREVIOUS  RUNS,  IT  WAS  FOUND  THAT  N1  IS  NOT  LESS 

0160  C  THAN  NSTAR/8.  THUS,  THE  INITIAL  VALUE  OF  N1  IS  SET  ACCORDINGLY. 

0161  C****************************************************** 

0162  C 

0163  120  DO  165  LL= INT ! NSTAR /8 ) , NSTAR 

0164  N1=LL 

0165  IF(N1. LT. C2 ) N1~R 1M1 

0166  C ************************ ****************************** 

0167  C  CHECK  B!N1.P0,  C2)  G. T.  1-ALPHA  CONSTRAINT 

0168  C ****************************************************** 

0169  125  CALL  PR0BS1 !N1» PO, C2,  BXLEC  > 

0170  IF!BXLEC.LT. ! 1. -ALPHA) ) GOTO  175 

0171  C****************************************************** 
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0172 
0173 
0174 
0175 
0176 
OJ  77 
0170 
0179 
0180 
0181 
0182 
0183 
0184 
0185 
0186 
0187 
0188 
0189 
0190 
0191 
0192 
0193 
0194 
0195 
0196 
0197 
0190 
0199 
0200 
0201 
0202 
0203 
0204 
0205 
0206 
0207 
0208 
0209 
0210 
0211 
0212 
0213 
0214 
0215 
0216 
0217 
0210 
0219 
0220 
0221 
0222 
0223 
0224 
0225 
0226 
0227 
0228 


C  CALCULATF  SECOND  SAMPLED 

C  ********************  ********  ************************ 

N2— NS7  AR  -Nl 

C ***********  *********  *********************************  * 

C  CHECK  THAI  DOUBLE  PROBABILITY  G. T.  1-ALPHA 

C ****************************** ************************ 

130  CALL  PR0BD1 (Nl, N2, PO, DPROB, KINDEX, R J ) 

IF ( DPRQB. GE. ( 1. -ALPHA) )GOTO  135 
IF-  (Cl.  EO.  C2M1  >  GOTO  175 
C1=C1+1 
GOTO  130 

C ************* ************** *************************** 

C  FIND  N2  VALUE  SATISFYING  AQQL  L. T.  AOQL* 

C ************************************ ****************** 

135  CAL  L  SRCH2(NNN,  Nl,  N2,  A0Q1  ) 

IF (A0Q1.LE. AOQL) GOTO  150 

C ****************************** ************************ 

C  N2  WILL  NOT  BE  LESS  THAN  Nl*9,  THE  INITIAL  LOWER  BOUND  ON  N2. 

C *************************** *************************** 

N2TEMP=N1*9 

C 

CALL  SRCH2(NNN,  Nl,  N2TEMP,  A0Q1  ) 

IF ( AOQ 1. GT. AOQL > GOTO  165 
C 

BL— N2 
BH“N1 *9 

140  N2TEMP=INT( (BL+BHJ/2) 

CALL  SRCH2LNNN,  Nl,  N2TEMP,  ADQ1 ) 

IF ( A0Q1 . LE. AOQL ) BH=N2TEMP 
IF  <  A0Q1 . LE. AOQL ) FA0Q=A0Q1 

C 

I F ( AOQ 1 .  GT. AOQL ) BL=N2TEMP 
IF< (BH-BL). EQ. 1. )GOTO  145 
GOTO  140 
145  N2=N2TEMP 

IFIAOQI. GT. AOQL )N2=BH 

C ****************************************************** 

C  CHECK  THAT  BINOMIAL  PROBABILITIES  ARE  G. T.  1-ALPHA 

C ****************** ************************************ 

CALL  PR0BD1 < Nl » N2, PO,  DPROB, KINDEX, R1 ) 

IF ( DPROB. LT . ( 1 . -ALPHA ) ) GOTO  165 
150  CALL  PROBSl(Nl, PO, Cl, BXLEC ) 

PA2=DPR0B-BXLEC 

C ****************************************************** 

C  COMPUTE  ATI 

C****************************************************** 

DAT I =N1 *DPR0B+N2*PA2+NNN* ( 1. -DPROB) 
C****************************************************** 

C  IF  THE  ATI  INCREASES,  CONTINUE  FOR  5  ADDITIONAL  INCREASING 

C  ITERATIONS.  THEN  INCREMENT  Cl  AND  CONTINUE. 

C ********************************* ********************* 

IF ( JJ. EQ. 4  >GQTQ  155 
IF(DATI. GE. DDATI > JJ=UJ+1 
IF (DAT I . GE. DDATI )GOTQ  165 
DDATI=DATI 
DDAOQ=FAOQ 
C 


36 


OAF  I 


0229 

Kl— C  1 

0230 

K2— C2 

0231 

K3=N1 

0232 

K4~N2 

0233 

C 

0234 

GQT  0  165 

0235 

C *#*•#■  (H*  ***■«■** **** *■*•****■**■**■*■* *•*■«•*#•*•** ***«•■«■ 

0236 

C 

MINIMUM  OF  CELL  (TCMIN)  FOUND 

0237 

c*«  *»***■**-****■*•*•*■**•#•****•***■•*  *************************** 

0238 

1  55 

TCMIN=DDAT 1 /NNN 

0239 

PDA! I=NNN 

0240 

C 

0241 

WRITE (5,  1 60 ) K 1  /  K2/  K3,  K4,  TCMIN 

0242 

WR  ITE  (  1 ,  160  )K1,  K2,  K3,  K4,  TCMIN 

0243 

160 

FORMAT!  1  OX,  12,  6X,  12,  5X,  14,  6X,  14,  7X,  F3.  6) 

0244 

C*<HHHf*******#*##*#*#**#*#******#*#***#it*************** 

0245 

c 

IF  MINIMUM  OF  COLUMN  IS  FOUND,  INCREASE  C2 

0246 

C  * ***«■•«•**•»**•»■** ********** ******* *•*■** **•****■*  **■*■#•  *#•»*•«■•** * 

0247 

IF ( TCMIN.  GE.  TTCMIN ) GOTO  170 

0248 

KK1=K1 

0249 

KK2=K2 

0250 

KK3=K3 

0251 

KK4=K4 

0252 

TTCMIN=TCMIN 

0253 

TT  AOQ=DDAOQ 

0254 

C 

0255 

IF  (C1.EQ.C2M1)  GOTO  175 

0256 

C1=C1+1 

0257 

C 

0258 

JJ=0 

0259 

N1=KK3 

0260 

c 

0261 

GOTO  125 

0262 

c 

0263 

165 

CONTINUE 

0264 

c****************************************************** 

0265 

c 

MINIMUM  OF  COLUMN  ( TMIN )  FOUND 

0266 

c****************************************************** 

0267 

170 

TMIN=TTCMIN 

0268 

TAOQ=TTAOQ 

0269 

C*****  ************************************************* 

0270 

c 

IF  MINIMUM  SAMPLING  PLAN  FOUND,  STOP 

0271 

c****************************************************** 

0272 

IFITMIN. GE.PLAN1G0T0  190 

0273 

c 

0274 

PLAN=TMIN 

0275 

KKK1=KK1 

0276 

KKK2=KK2 

0277 

KKK3=KK3 

0278 

KKK4=KK4 

0279 

175 

C2=C2+ 1 

0280 

c*********************»******************************** 

0281 

c 

NEW  BOUNDS  ON  SAMPLING  PLAN  CALCULATED  FOR  NEW  VALUE 

0282 

C  *•#■******•**«•**  **#***#####*******#**# fl-**#******##***#*** 

0283 

180 

CALL  SEARCH ( NNN,  C2,  NSTAR,  AOQ) 

0284 

IF  < AQQ. LE. AOGL ) GOTO  185 

0285 

IF ( NSTAR. GT.  NNN) GOTO  175 
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0286 

0287 

0280 

0289 

0290 

0291 

0292 

0293 

0294 

0295 

0296 

0297 

0298 

0299 

0300 

0301 

0302 

0303 

0304 

0305 

0306 

0307 


C 

c 


c 

c 


NS  f  AR'NST  AR+-1 
GOTO  180 


185  CALL  PR0BS1  <  NSf  AR  .  PO  C2  BXI.EC) 
GOTO  110 


J  90 
1  95 


?00 


205 


210 

215 


WRITE (5. 195) 

FORMAT ( / / 1  OX ,  'SAMPLING  PLAN  MINIMUMS') 

WRITE! 5, 200 ) KKK 1  ,  KKK2 

WRITE! 1- 200 ) KKK 1,  KKK2 

FORMAT  !  /  1  OX >  'Cl«  ',  12,  2X,  'C2=-  ',12) 

WR  I TE  !  5,  205  )  KKK3,  KKK4 
WRITE! 1, 205 ) KKK3,  KKK4 
FORMAT!  10X,  'Nl-  13,  2X,  'N2-  ',  13) 

WR ITE( 5* 210) PLAN 

WRITE!  1, 210)PLAN  ,  N 

FORMAT ! 10X ,  'MINIMUM  AFI=',T8.6) 

RETURN 

END 
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0001 
0002 
0003 
0004 
0005 
0000 
0007 
OOOO 
0009 
0010 
OOl  1 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 


SUDROUT INE  SEARCH ( NNN,  C .  NS,  AOQ ) 

C *******  ##*#*»***•«****■***■* **##*****-tl--tt ***************-*<  «-* 
C  SEARCH  ro  FIND  VALUE  OF  PS  TAR  USING  GOLDEN 

C  SECTION  METHOD.  INITIAL  LIMITS  O'  O  AND  1 

C  tt*#*#*#*****'*#*'**-*-*-*#-#-#*#-*****  ,:  «■«■# 

INTEGER  C 

DOODLE  PR  I  C,  J  SION  SUMLOO 
CQMMQN/DI  K7/SUML0G ( 4000) 

COMMON /ni  KD/N 

c: 

A 1-0. DO 
A2--  1 .  DO 
11  .  D-3 

R- 5. D- 1 M DSORT ( 5. DO ) - 1 . DO ) 

H -A2-A1 

PLFT  =A  1  +  ( Rc-K  )  *H 
PRT -A  1  +  <  R*H ) 

C 

C 

CALL  PR0DS1 (NS, PLFT, C, BXLEC) 

AT  I = ( NS*BXLEC ) +NNN* ( 1. DO-BXLEC ) 

AF I  1 =AT I /NNN 
A0G1=PLFT*( 1. DO— AF I  1 > 

CALL  PRODS 1 (NS, PRT,  C.  BXLEC  > 

AT  1=  (  NSK-BXL.EC  )  +NNN*  (  1.  DO-BXLEC) 

AF I 2=AT I /NNN 
A0G2=PRT*  (  1 .  D0  -AFI2  ) 

GOTO  110 
C 
C 

lOO  CALL  PR0DS1 (NS, PLFT, C, BXLEC ) 

AT I = <  NS*B  XLE  C ) +NNN*  < 1 . DO-B  XLFC ) 

AF 1 1 =ATI /NNN 
A0Q1=PLFT  M 1. DO-AFI 1 ) 

GO  TO  110 
C 

105  CALL  PR0BS1 (NS, PRT, C, BXLEC) 

ATI  =  (NS*BXLEC  > +NNN* ( 1. DO-BXLEC) 

AF I 2=AT I /NNN 
AOQ2=PRT*( 1. DO  -AFI2) 

C 

HO  I F  ( AOQ  1.  l-T.  A0Q2)  GOTO  115 
A2— PRT 
H=PRT-A1 

IF ( ABS ( PRT  PLFT ). LE. T ) GOTO  120 
C 

PRT=PLFT 
PLFT=A1+(R*R)*H 
AOG2-AOQ1 
GO  TO  100 
115  A 1 =PLFT 
H-A2— PLFT 
C 

IF ( ABS ( PRT  PLFT). LE. T ) GOTO  120 

PLFT=-PRT 

PRT=A1+R*H 

A0Q1=A0Q2 
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0053 

0057 

0060 

0061 

0062 

0063 

006*1 

0065 

0066 


GO  TO  105 

C 

C 

120  PS~  ( PL  F  1  t-PRT  )  /2.  DO 

CALI.  PRUBS1 (NS, PS, C, BXLEC > 

AF  I  —  <  (  NS*BXLEC  )  +NNN  *  (  1.  DO  BXt  EC  )  )/NNN 
AOG-PS* < 1. DO  AF 1  ) 

RETURN 

END 
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0001 

0002 

0003 

000 A 
0000 
0006 
0007 
oooo 

0007 
0010 
00  t  1 
0012 
oo  i  ;i 

00)  4 
OO  .1  5 
00 1 6 
00  1 7 
OOlIt 
OO  1  7 
0020 
002  1 
0022 
0022 
0024 
0020 
0026 
0027 
0028 
0027 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0037 
0040 
0041 
0042 
004  3 
0044 
0045 
0046 
0047 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
005/ 


SUB  ROUT  I  N!_:  0HCH2(NNN,  N1  -  N2,  AUO  > 

C 

INTEGER  Cl,  02.  IM 
COMMON/BI  K6/C1,  (  2 

C 

K INDEX  I 
A 1  O. 

A2  1  . 

I  0001 

R  - .  5MDSGRI  (5.  DO;  1.  ' 

IT- AO  A  1 

Pl_l  r  -  A1  HIB-H  )  (HI 
PRT  -A  1  MR  (HI) 

0 

C 

CALI .  PRUBH1  <Ni,  PLFT  .  Cl,  PAl  ) 

CALI.  PR0BD1  <  N  1 .  M2.  PL.FT  .  DPRDB,  K  INDEX.  R  I  ) 
PA2--DPR0B  PAl 

ATI  —  DP ROB *N  1  -*  P  A2*N2+NNN)  (  1  -DPROB) 

API  1 =AT I /NNN 

AGO 1  - PLFT (  1  .  -API  1  ) 

CALL  PR0BS1 <N1, PRT, Cl, PAl ) 

CALL  PR0DD1  <N1  ,  N2.  PRT,  DPROB,  KINDEX,  R1  ) 
PA2=DPR0B  -PAl 

ATI=DPR0B*Nl+PA2*N2+NNN*-(  1,  -  DPRC  B  ) 

AFI2-ATI/NNN 

AOQ2—PR7 * ( 1 .  -APJ2) 

GOTO  110 

C 

c 

lOO  CALL  PR0BSUN1,  PLFT,  Cl,  PAl  > 

CALL  PRQBD1 <  Nl , N2, PLFT,  DPROB, KINDEX,  R1  ) 
PA2=DPR0B  -PAl 

ATI=DPR8B*N1+PA2*N2+NNN*( 1  - -DPROB ) 

AFI 1=ATI /NNN 
AOQl^PLFT * ( 1 .  -AFI 1 ) 

GOTO  110 

105  CALL  PR0BS1 <N1 , PRT, Cl, PA)  ) 

CALL  PR0BD1 (Nl,  N2,  PRT, DPROB, K INDEX,  RJ  ) 
PA2=DPR0D  PAl 

ATI=DPR0B*IM1+PA2*N2+NNN*< 1. -DPROB) 

AFI 2— ATI /NNN 
AOQ2=PRT*<  1.  -Af  12) 

C 

C 

110  IFIAOQI. L T. A0Q2)G0T0  115 

C 

c 

A2:  tPRT 
H-PRT-A1 

IF  (  ABS  (  PR  1  -PLFT  ).  LE.  T)  GOTO  120 
PRT  =PL FT 
PLFT=A1  MRc-R  )  *H 
A0G2- AQG 1 
GOTO  100 
115  A1 -PL.FT 
C 
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rr^ 

c’c 

oc 

c 

il-AL'-PLI  1 

0060 

lFLABSIPPI  PIT  f  ) . Lh  .  r  ’>  GO  1  0  (30 

006  1 

pi. r  r  •  pr  i 

006;' 

PR  t  --A1+R  *•  H 

0063 

AGO 1 - AGO? 

0064 

0060 

(. 

GOT  (.1  105 

0066 

0067 

c 

130  PS  ( PL.F 1 +PR 1)/?. 

006(3 

CALL  PROBS1 <N1,  PS,  Cl  PA1  ) 

0069 

CALL  PR0BD1 ( N1 ,  N2  PS,  DPRGB ,  K I MDEX ,  R 1  ) 

oo/o 

HA2DPR0B  PAL 

0071 

AF  I  —  (  DPROB  M-Nl  +PA2*N2->NNN*  (  1.  -  DP RGB  )  )  /NMM 

00/? 

A0Q=PS*(1.  AFT) 

0073 

RETURN 

0074 

PND 
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OOO ) 

0002 

0003 

(>00«1 

OOO') 

000* 

0007 

ooom 

OOO'/ 
00  JO 
00.1  J 
001  2 
OO  J  3 
00 1  4 
00  1  0 
00  1  * 
0017 

oo  i  e 

0019 

0020 

0021 

0022 

0023 

0024 

0020 

002* 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

003* 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


SUHROUT  INI  :  ’RUH3  1  (NN,  P,  0  .  B  X  L  E  (  ) 

(>«#***(*  ***•»#*  t 

C  THIS  SUB  ROUT  I  Nf  COMPUTES  CUMULATIVE  BINOMJAt 

C  PR01IADIL1  I  It  S 

(:»«  »  !>#**»«#<**««*»* 
INTEGER  C 

DOUBLE  PRECISION  SUMLUG 

C 

COMMON/BI  K7/SUMI  UG  (  4000) 
coMMON/m  k;j/n 

( 

0  J  .  -P 

(:»»****»»**##««»*♦**♦**«****»#*#*«*#♦**»«««*»*♦*«**#*#< 

C  BINOMIAL  PROB.  WHEN  C-0 

CfHHHHH  IHHt****  ►#»*«»*»»*»*»**<■**♦**«««*  »****■«  ****•»••«■■«•* «-* 

CSUMS -0«  «NN 

IE  <C, EG. O)  GOTO  45 

C  *»*»»«♦■«■*♦♦*#  fr#***#*******#**##-*  »*■!/«-*«■**  *■»*«■#  *■«■* 

C  AVOID  RECOMPUTING  SUMLOG(I)'S  ALREADY  IN  MEMORY 

C#*#»ft*##*#***  «■**  **•»■•»■■»■*•»■*•»■•»■•»■*■»■*■«•■«■  ♦■ft******************** 

IE  (  N-  NN  )  10,25,25 
JO  M-N+i 

C  ********•***■*■»■■*■*■*  a-******-#-  *******  ***■«•  «■*•«•  *#*  **  *  *■**■*  **■•*■**  * 

C  COMPUTE  N  SUMLOGS-EQU I VALENT  TO  N-FACTOR I AL 

C#*»lH>**»**iHHt#*»#**#*******«***#***»*<H(*#*»**#**#*)Ht** 

IF  (M. GT. 1 )  GOTO  15 
SUMI.OG  (  1  >-0. 

IF ( NN. LE. 1 )  GOTO  25 
M-2 

15  DO  20  I=N, NN 

SUMLOG ( I >  =DLOG 1 O ( DFLDAT ( I ) ) +SUMLOG ( I - 1 ) 

20  CONTINUE 

CiUHt*##*##*#******#**###**##*#*#**********#****##***#** 

C  COMPUTE  C  SUMS-EGU I VALENT  TO  SSUM  OF  PROB. COMP  IN. 

C  I.E.  CUMULATIVE  BINOMIAL  DISTRIBUTION  COMPUTATION 

C#*«***##*****fr***#*****###******HHHHt******#****#»**### 

25  IF (NN. GT. N)  N=NN 

C  ****■«■*■•*■***■*•*'**«•#*•**■***■*■**•*■#■#***■»■*■*•&**•*•*****-****#  ****** 

C  DETERMINE  BEST  NUMBER  HANDLING  LOOP 

C*****#####****-#*###*##*-**#*#******#*****##*##****#**** 

IF  (NN. GT.300)  GOTO  35 
DO  30  K=l,  C 

CSUMS~ 10. **< SUMLOG ( NN ) -SUMLOG ( NN-K ) -SUMLOG ( K ) ) 

1  *P**K*Q**(NN-K> vCSUMS 

30  CONTINUE 
GOTO  45 

C  ***•»•«•#******  ******  ******  *******•«•***■•**■**■#***■***«•*##*#*■* 

C  LOOP  FOR  LARGE  EXPONENTS 

C##*#**###'*'**#*#*##*#**###**###*#*#*#*#*****#**#*****#* 

35  DO  40  K=1 ,  C 

CSUMS- 10. ** ( SUMLOG ( NN ) -SUMLUG ( NN-K ) -SUMLOG ( K ) 

1  +K*DLOG 1 0  <  DBLE  <  P  >  >  +  <  NN-K ) *DLOG 1 0  <  DBLE ( Q ) ) ) +CSUMS 

40  CONTINUE 

C 

45  BXLEC  ^  CSUMS 
RETURN 
FND 
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SUBROUTINE  PR0BD1 <N1,  N2,  P.  DPROB/  K,  HI  > 

C  ********  *#**###*##***#*****#*■»****«■***##•»#**********•«■* 

C  THIS  SUBROUTINE  COMPUTES  DOUBLE  PROBABILITIES  FOR 

C  COMPUTING  SECOND  SAMPLE  NUMBER  OF  DOUBLE  SAMPLING  NUMBER 

C  ******************************************************** 
CQMMQN/BLK6/C 1 ,  C2 
INTEGER  Cl . C2- R1 
C 

IF  ( K. EQ. 1 )  CALL  PR0BS1 (Nl. P, Cl. BXLEC ) 

IF (K. EQ. 2 )  CALL  PR0BS2 ( Nl , P , C 1 . BXLEC ) 

DPROB=B  XL.EC 
TEMP-BXLEC 
NTEMP=C  1  1 

KTEMP-R1-1 

DO  10  IX-NTEMP, KTEMP 
I  =  IX 
J=C2-I 

IF(K.EQ. 10  CALL  PROBS1 ( Nl , P, I , BXLEC ) 

IF ( K. EQ. 2 0  CALL  PR0BS2 < Nl . P, I , BXLEC > 

PROB 1=BXLEC— TEMP 
TEMP=BXLEC 

IF ( K. EQ. 1 0  CALL  PROBS1  <  N2,  P.  J,  BXLEC  > 

IF (K. EQ. 2)  CALL  PR0BS2  ( N2,  P.  U,  BXLEC  0 
DPRQB=DPROB+  <  PROB 1 *BXLEC > 

10  CONTINUE 
C 

RETURN 

END 
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